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ABSTRACT 


This document gives calculations from the best avail- 
able (summer 1964) geodetic sources showing that the Syn- 
com II - Brazil second-order earth gravity experiment 
reported earlier by this author probably contained a small 
but significant bias due to higher order earth gravity. This 
bias is estimated by "extrapolation" from these sources. 

Bounds on the true ellipticity of the earth’s equator are 
derived by combining the best estimate of this bias with the 
actual Syncom II drift measurements. These bounds, though 
somewhat wider than measured by Syncom n (over Brazil) 
alone, should be absolute and the target for all future in- 
vestigations into longitude-dependent earth gravity. 
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ON THE PROBABLE INFLUENCE OF HIGHER ORDER 
EARTH GRAVITY IN THE DETERMINATION 
OF THE EQUATORIAL ELLIPTICITY OF THE EARTH 
FROM THE DRIFT OF SYNCOM II OVER BRAZIL 

by 

C. A. Wagner 

Goddard Space Flight Center 

INTRODUCTION 

Purpose 

This report makes use of lower altitude gravity data to evaluate the likely bias due to higher 
order earth effects in the Syncom n - Brazil second-order gravity-drift experiment. 

The Syncom II - Brazil drift, previously analyzed for second-order longitude gravity only (i.e., 
the ellipticity of the earth’s equator) sampled a limited portion of that longitude field. From this 
limited sampling, no direct evidence of higher order earth gravity effects at 24-hour satellite 
altitudes was available. Nevertheless, preliminary calculations using the evidence of a number of 
"full field" geoids derived from lower altitude gravity observations appeared to indicate that the 
assumption (in the Syncom II gravity experiment) of small higher order gravity effects at 24-hour 
altitudes was justified. The reason for this study is to present detailed calculations "extrapolating" 
the evidence from lower altitude gravity surveys in support or contradiction of this judgment * 

If the effects of higher order earth gravity at 24-hour altitudes are large and not well deter- 
mined from the lower altitude evidence, then the isolated accelerations measured by Syncom n 
drift over Brazil cannot fairly be extrapolated to provide comprehensive information about the 
likely 24-hour longitude gravity accelerations at other longitudes. However, if, as assumed in the 
previous study, Syncom II drift accelerations over Brazil were likely to be fairly representative 
of the second-order longitude gravity field only, then these isolated but well- measured accelera- 
tions tell us significant information about the probable drift accelerations on 24-hour satellites at 
all longitudes. 

Background on the External Gravity Field of the Earth 

For over 100 years scientists have represented the earth’s actual external gravity field as an 
infinite series of spherical harmonic components (Appendix A, Table A-l), each representing a 


•Subsequent 24 hour satellite data has confirmed the minor influence of higher order gravity on the Syncom -Brazil measurements, as 
predicted in this report. 
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small distortion of the dominant spherical field first calculated by Isaac Newton in the 17th century 
(Reference 1). It can be appreciated that every measurement of this field at a single point in geo- 
graphical space involves all the terms in this series. Thus, theoretically, one can never deter- 
mine each term of the series with complete certainty unless one has an infinite number of exact 
point measurements to deal with. However, as a practical matter, the earth proves to be so nearly 
spherically homogeneous that the first-order term of the gravity series has changed but little sinc< 
Newton's time in the 17th century. Nor is there any reason to suppose it will drastically alter in 
the future when even more gravity data will become available. The same hopeful remarks can be 
made for the dominant second-order gravity component that nearly perfectly describes the polar 
flattening of the earth due to the relatively strong centrifugal action of its rotation. For evidence 
of this the reader is referred to the many recent studies of the gravity field in the literature which 
report values for this "oblate” term, particularly those which contain sufficient data to allow for a 
simultaneous determination of a large number of higher order "oblate" terms (Reference 2, for 
example). A more complete discussion of this problem with respect to the oblate and zonal terms 
of the gravity field is found in Reference 3. With respect to the very small longitude components 
of the field, the situation is far more difficult to resolve. This is so, essentially because the earth 
does not seem to be able to support mechanically large mass inhomogeneities in the upper mantle 
and crust. One can avoid the problem of finding approximations to the true longitude field valid 
everywhere by reporting effectively the point force measurements themselves and leaving the stud 
of what may appear to be longitude gravity "noise" to some future date. Nevertheless, there are 
theoretical expectations (associated with the deep earth convection currents that are supposed to 
be the cause of the external magnetic field) for a dominant longitude gravity field distortion of low 
order. But the search for a well-defined, lowest (second) order longitude field, represented by an 
earth with an elliptical equator, had been disappointing and inconclusive until the advent of Syncom 
n, the first operational 24-hour satellite (Reference 4). The great height of this satellite makes it 
sensitive mainly to gravity terms of a low order associated with mass distortions distributed over 
wide areas of the earth's surface. 

Method of Analysis 

Anomalous long term accelerations of Syncom n were heretofore assumed by this author 
(Reference 4) to be caused by the second-order longitude gravity field of the earth only. Since 
essentially only two independent point measurements of the field were obtained from the limited 
extent of drift data in Reference 4, no information on higher order distorting effects can be ob- 
tained without recourse to gravity data from other sources. 

This study selects eight recent and largely independently derived "full" fourth-order geoids 
as the basis for assessing the likely influence of higher than second-order earth gravity on the 
Syncom II - Brazil drift accelerations. Over each geoid the full field drift accelerations on the 
two Syncom n - Brazil orbits are calculated. As in the actual Syncom n - Brazil experiment, 
these calculated full field accelerations for each geoid are then assumed to arise from a simple 
second-order field only. This "two point" second-order longitude field (essentially giving the 
magnitude and longitude orientation of the earth's elliptical equator) is then calculated. Direct 
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cQmparison is then made of this "two point ,T second-order field with the "actual" second-order 
field for each geoid. This comparison gives the bias implicit in the T *two point" field measured 
under Syncom n - Brazil orbit conditions for each fourth-order geoid. It is assumed that earth 
gravity effects higher than the fourth order have negligible effect on the 24-hour satellite. 

A SYNCOM II ORBIT GRAVITY EXPERIMENT 
SIMULATED IN A FOURTH-ORDER FIELD 

Imagine the Syncom n longitude sampling over Brazil (Reference 4) to have taken place in the 
earth gravity fields of Appendix A. (See also Figure 1.) From Reference 5 it is shown that the 
orbit-averaged longitude acceleration of the 24-hour satellite figure eight ground track (or ascend- 
ing equator crossing of the satellite) is given as 


X = -12 TT : 
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Equation 1, with only the j 22 term, was first derived in Reference 4. The right-hand side is es- 
sentially the orbit- averaged tangential (energy-changing) force on the 24-hour inclined-orbit satel- 
lite. The left-hand side is the longitude acceleration caused by the continual period change effected 
by this force. 


3 


ORBIT AVERAGED EARTH GRAVITY LONGITUDE 
ACCELERATION ON A 24- HOUR 32.9° INCLINED 
SATELLITE (units: -12ir2 (R 0 /a,) 2 X 10' 5 rad/(sid day) 2 ) 
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Figure 1— Full-field orbit-averaged longitude accelerations on Syncom II as predicted 
by six recent geoids and as actually measured in a Syncom II gravity experiment. 

The inclination of Syncom n during the excursions over Brazil reported in Reference 4 was 
i = 32.9°. With this value in Equations 2-6 the "inclination factors" for Syncom n in Aug. - Feb. 
1963-64 are: 

F( i ) 22 = 0.845 

F( i ) 31 =0.267 

F( i ) 33 = 0.775 

F( i) 42 = 0.046 

F(i) 44 = 0.710 . (7) 
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From Reference 4, (R 0 /a s ) = (6378/42166) = 0.1513 and (R 0/ /a s ) 2 = 0.0229. Equation 1 then 
gives 

X = (5.06J 22 sin 2(\-\ 22 ) -0.0605 J 31 sin(\-\ 31 ) + 5. 26 J 33 sin 3^ - \ 33 ) 

-0.0158 J 42 sin 2(a. - A. 42 ) + 6. 83 sin 4(A - X- 44 )} ' (») 

which is in units of - 12 tt 2 (R 0 /a s ) 2 rad/day 2 . Figure 1 gives a graph of Equation 8 for six recent 
geoids. For the first "synchronous” longitude location of Syncom EE at k = - 54.8°, Equation 8 gives 

X_ S4 g . = ^5.06 J 22 sin 2(-54.8° -\ 22 ) -0.0605 J 31 sin (-54. 8° - A- 31 ) + 5.26J 33 sin 3(-54. 8° - \ 33 ) 

-0.0158 J 42 sin 2(-54.8° -\ 42 ) + 6.83J 44 sin4(-54-8° -\ 44 )j * (9a) 

For the second "synchronous" longitude location of Syncom n at A. = - 59.2° Equation 8 is 

A.. r = { 5 . 06 J 22 sin 2(-59. 2° -A 22 ) -0.0605 J 31 sin (-59. 2° - \ 31 ) +5.26J 33 sin3(-59.2° - 

-0.0158 J 42 sin 2(-59. 2° -X 42 ) + 6.83J 44 sin4(~59.2° -\ 44 )} * (9b) 

In the actual gravity experiment on Syncom n in Reference 4 it was assumed that Equations 9a and 
9b could be truncated at the first right-hand terms; or that longitude gravity terms higher than 
second-order had negligible net effect on the longitude acceleration. We can test this assumption 
by performing the identical experiment with the "fourth-order field" accelerations of Equations 9a 
and 9b (as assessed by the geoids of Appendix A) under the same assumption. 

For each geoid a j 22 and k 22 is listed in Table 1. These values are assumed to be best for a 
certain set of gravity data taken over a large zone of geographic space. If one assumes in each 
case that the high-order geoid is perfect, then apparent J 22 and k 22 values can be computed from 
a limited two- point measurement of the field of that geoid as in the actual Syncom n experiment 
over Brazil. Comparing the two sets of coefficients gives a direct measure of the distortion in j 22 
and k 22 as it is computed in the Syncom II - Brazil experiment which would be expected if that 
particular geoid were perfect (see DISCUSSION). 

For each geoid, a J 22 apparent (j 22 ) and a k 22 apparent (\ 22 ) may be computed from the fourth- 
order field accelerations of Equations 9a and 9b under the simplifying assumption of Reference 4. 
Thus X (-59.2°), as it is calculated from Equation 9a, is the fourth-order field acceleration which 
would have been measured by Syncom n at -59.2° for the geoid whose gravity constants are J nm and 
k na . Under the simplifying assumption of the actual Syncom n - Brazil Experiment (Reference 4), 
this higher order acceleration (actually measured) was assigned to only the second order gravity 
field, giving apparent earth constants j 22 and k 22 . 
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Table 1 


Higher Order Gravity Effects on Equatorial Ellipticity as Seen 
by a Syncom II - Orbit at Two Longitudes, for Eight High-Order Geoids. 





Full Field Acceleration on 
Syncom II, Orbit Averaged 
(see Figure 1) 

(units of -127T 2 (Ro/ a s ) 2 x 10 “ 6 rad/(sid day) 2 ) 

J 22 (reported) 


Geo id 

J 22 (reported) 

k^ (reported) 
(degrees) 

^22 ” J 22 (apparent)*** 
(\ = - 54.8° and 

- k 22 (reported) 
-^22 (apparent)*** 




\ = - 54.8° 

k = - 59.2° 

k - - 59.2°) 


Kaula-Comb. (1964) 

-1.51 xlO' 6 

-15.5° 

6.793 

6.815 

1.12 

-2.9° 

Izsak (1964 

-1.00 

-17.0° 

3.876 

4.099 

1.19 

+5.0° 

Kaula (1964) 

-1.77 

-18.2° 

7.962 

8.194 

1.085 

i 

© 

CO 

o 

Guier (1963) 

-1.80 

-10.4° 

8.115 

7.944 

1.12 

-2.3° 

Uotila (1962) 

-1.52 

-36.5° 

3.130 

3.682 

1.55 

-1.2° 

Kozai (1962) 

-1.2 

-26.4° 

4.287 

4.809 

1.07 

+3.9° 

Zhongolovitch 

(1957) 

-5.95 

-7.7° 

28.684 

27.436 

1.025 

o 

00 

CO 

1 

Jeffreys (1942) 

-4.1 

0.0 

18.184 

16.888 

1.06 

1 

© 

bb 

o 

Actually Measured 
in Syncom II - 
Brazil Experiment 



8.1* 

8.45** 

Mean = 1.095 
(excluding 
Uotila (1962)) 

Mean = -0.4° 
sample root 
mean square 
deviation is 2.9° 


♦From Reference 4, X(-54.8°) = - 2.20 x 10~ 5 rad/(sid day) 2 
♦♦From Reference 4, M-59.2°) = - 2.29 x 10" s rad/(sid day) 2 


♦♦♦Calculated values from the full field accelerations of the k - ~ 54.8° and k - ~ 59*2° columns, assuming these accelerations 
arise from a J 22 field only. 


Duplicating the experimental procedure for each of the geoids in Appendix A, we set Equation 9a 
to 


^- 59 . 2 - = 5.06J 2 ' 2 sin2(-59.2°-\ 2 , 2 ) 


5. 06 sin (-118. 4°) [j 2 ' 2 cos ] + 5. 06 cos (-118.4°) [- J 2 ' 2 sin 2\,' 2 ] 


(10a) 
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Similarly, we simulate the actual Syncom II - Brazil experiment at \ = - 54.8° by setting Equa- 
tion 9b to 


^-54.8- = 5.06 J 22 sin 2(-54.8°-\ 2 ' 2 ) 

i 

= 5.06 sin (-109.6°) [j 2 ' 2 cos 2^ 2 ’ 2 ] + 5. 06cos (-109.6°) [~J 2 ' 2 sin 2^j’ 2 ] (10b) 

I 

Solving Equations 10a and 10b for (j 22 cos 2^ 22 ) and (- J 22 sin 2A 2 ' 2 ) in terms of \(-54.8°) and 
X-(-59.2°) for a fourth-order field will give 


and 


J 22 COS 2A-22 = 


A.(~59. 2°) 5.06cos (-118.4°) 

X(-54.8°) 5.06 cos (-109 6°) 


5.06 sin (-118.4°) 5. 06 cos (-118. 4°) 

5.06 sin (-109.6°) 5.06 cos (-109.6°) 


{X(-59.2°) cos (-109.6°)-l(-54.8°) cos (-118.4°)} 
[5.06 sin (-8.8°)] 


(Ha) 


- J ' sin 2K' 


|5. 06 sin (-118.4°) M~59. 2°) 

' 

5.06 sin (-109.6°) *-(-54.8°)| 


(5.06) 2 [sin (-8.8°)] 


{£(-54- 8°) sin(-118.4°)-X(-59.2°) sin (-109. 6°)} 

5.06 [sin (-8.8°)] 


(lib) 


When Equation 11a is divided by Equation lib, \ 22 (apparent) may be calculated as in the actual 
Syncom n - Brazil experiment, from 


K 


22 


1 

2 tan 


~ ['((-54. 8°) sin (-118.4°) -X(~59.2°) sin (-109.6°)] 
A.(-59.2°) cos (-109.6°) -\(-54.8°) cos (-118.4°) 


(12a) 
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Similarly, when Equations 11a and lib are squared, J 22 (apparent) may be calculated, as in the 
actual Syncom II experiment, from 


J 22 


- 1(T 6 

5.06 sin (8.8°) 


— 54. 8°) sin (-118.4°) -X(-59.2°) sin (-109.6°)] : 


+ [X(-59. 2°) cos (-109.6°) -K(-54.8°) cos (-118.4 0 )] 2 j 1/2 • (12b) 


When the following values are used: 


sin (-8.8°) = -*0.15299 

sin (-118. 4°) = -0.87965 

cos (-118.4°) = -0.47562 

sin (-109. 6°) = -0.94206 

cos (-109.6°) = -0.33545 , 

Equations 12a and 12b become 


\ 




0. 87965 '^(“54. 8°) + 0.94206 ^(”59. 2° 
0.33545 \(-59.2°)+0.47562\(-54. 8° 



(13a) 


and 


j; 


22 


1.292 x 10" 6 


[[-0. 87965 X(-54. 8) -0. 94206 X(~59 . 2)] 2 
| + [-0. 33545 X( -59. 2) + 0. 47562 X.( ”54. 8)] : 


(13b) 


For each geoid in Appendix A one can calculate a ratio R 22 = j 22 /j 2 ' 2 from Equation 13b, giv- 
ing the gross higher order effect on the magnitude of equatorial ellipticity which would have 
been seen by Syncom II if that geoid were correct. The cumulative divergence of these geoids 
from "reality" at the two Syncom n longitudes can best be seen in Figure 1. The cumulative 
gravity forces of the Syncom II - Brazil geoid (j 22 = -1.7 x 10“ 6 , \ 22 = -19°) at those longitudes 
are exact to within about 3%, since the full field forces were those actually sensed in the ex- 
periment of Reference 4. 
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Similarly, for each geoid in Appendix A one can calculate a major axis change, &k 22 = k 22 (for 
that geoid) - ^ 2 ' 2 (from Equation 13a), giving the gross higher order effect on the location of the 
equatorial major axis which would have been seen by Syncom n if that geoid were correct. 

In Appendix B this calculation, which is just a simulation of the Syncom n experiment in eight 
recent gravity fields reporting individual higher order tesseral coefficients, is carried through 
for the geoid of Kaula-Combined (1964). The results for this higher order geoid and the others of 
Appendix A are presented in Table 1. 

PROBABLE RANGE OF HIGHER ORDER DISTORTIONS 
TO EQUATORIAL ELLIPTICITY AS SEEN BY SYNCOM II OVER BRAZIL 

From this simulated study, reasonable bounds on the bias to be applied to the actual Syncom II 
determination of earth ellipticity appear to be 

R 22 = 1.1 ± .1 

and 

M 22 = - ( 1 ± 3)° 

The "mean" values of the bias factors in Equations 14 are, together, close to the factors computed 
for Kaula (1964) and Guier (1963) which give the full field accelerations for Syncom II closest to 
those actually observed (Table 1 and Figure 1). 

The Syncom n reduction of the observed full field accelerations, assuming only the second 
order effect, gave (Reference 4) 

J 22 (apparent) = - (1.70 +.05) x 10~ 6 

and 

\ 22 (apparent) = -(19±6)° • 

When the bias factors of Equations 14 are applied to Equations 15, the j 22 and \ 22 values from the 
Syncom n - Brazil experiment corrected for higher order earth effects become 

J 22 = - (1.70 ±0.05)(1.1 ±0.1)x 10' 6 = - (1.87 +0.18)x 10~ 6 (16a) 

^22 = - [(19±6)°-(1±3)°] = -(20 ±6.5)°. (16b) 
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It is to be especially noted that though the full field forces (or accelerations) at the two Syn- 
com II longitudes are widely different for the eight geoids, the distorting effects on the simulated * 
j 22 reductions, assuming a second-order field only, are small and consistent. For all geoids 
tested, the result of a minimum sampling of two longitudes of the field at 54.8° West and 59.2° West 
was to yield values of j 2 ' 2 or j 22 (apparent) up to 20% less than the reported J 22 value for that 
geoid. The anomalous result from Uotila (1962) is due to an overemphasis of the third-order 
harmonics in that geoid which is not shown by either the recent satellite geoids or the older gravi- 
metric geoids. The distortions in ^ 22 are all less than 5.0° and appear to have a westward bias of 
the order of 1 °. 

DISCUSSION 

The calculation of probable higher order distortion of true equatorial ellipticity as seen by 
Syncom n over Brazil given in the first two sections of this paper, did not utilize directly the con- 
straints of the observed accelerations of the satellite. It was felt that, except in two cases, the 
observed accelerations were so far from those predicted by the lower altitude geoids that to at- 
tempt to force fit observed accelerations with those predicted from the lower altitude data would 
be no more meaningful than the straightforward simulation actually calculated. 

These simulations at least have the virtue of assessing precisely and directly the distortions 
of "reality" we are interested in finding. Of course the "real earth" in each simulation is only an 
estimate, good or poor as the case may be. Yet there is encouraging consistency in the "internal" 
distortion results calculated for such widely differing geoids. By way of comparison, Kaula (Ref- 
erence 7) has performed a calculation for the second-order gravity field constrained by the Syn- 
com n accelerations directly. In this calculation all higher order effects were accumulated as a 
single third-order effect whose relationship with the second- order one was established by an 
averaging of the results of recent geoids. The conclusion was similar to the one obtained here. 
Kaula obtained revised Syncom n values of 

J 22 - - 1.87 x 10" 6 


and 


\ 22 = - 22.3° , 

giving a small strengthening and westward shifting of equatorial ellipticity over that found by ignor- 
ing higher order effects altogether. 

The method of estimation employed here is felt to be more comprehensive, to include a wider 
range of effects, and to be free from any preassignment of the relationship between them. 

The geoids chosen are felt to be fairly representative, without duplication, of the great amount 
and variety of longitude gravity data that has been accumulated over the years (Table A-l). The 
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geoid selection in each case but one was made on the basis of its being the latest and presumably 
the most comprehensive higher than second-order geoid by that author. The geoid Kaula- Combined 
(1964) was chosen as an example of a weighted estimate based on recent satellite-camera, satellite- 
Doppler and surface gravimeter geoids. Of the eight geoids, three deal with satellite-camera data 
(Izsak (1964), Kaula (1964) and Kozai (1962)) reduced from the arcs of up to ten medium altitude 
and medium inclination satellite orbits. Izsak (1964) and Kaula (1964) include reductions from a 
number of high inclination orbits as well. Guier (1963) is an example of a recent geoid reduced 
from saiellite-Doppler data cm a number of "TRANSIT” satellites. 

The three surface gravimeter geoids are the most recent (Uotila (1962)); a recent Russian 
geoid (Zhongolovitch (1957)) which was chosen for comparison because both the data and the 
statistical reduction techniques behind it are significantly different from those used as the basis 
for the other two gravimetric geoids; and the oldest geoid reporting terms higher than second 
order, which is by a highly regarded authority in the field (Jeffreys (1942)). 

In any case, it can be confidently expected that further precise Syncom II drift data at other 
longitudes will tell geodesists much more about the true extent of higher order longitude gravity 
than the rough bounds which have been established here by dealing with imprecise geoids derived 
from lower altitude data. The important fact to be gleaned from this study is that the bias in the 
simple second-order determination from Syncom II observations due to higher order effects is, 
from all indications, small. 

In Equation 16b the estimated total deviation in k 22 is taken as the square root of the sums of 
the squares of the component deviations, as is customary in the theory of errors for the sum of 
two random variables. In Equation 16a the estimated total deviation in j 22 is taken as the square 
root of the sum of the cross products of squared means and deviations (Reference 6). Rounding 
off these results upward to maintain a conservative estimate, we estimate the realistic bounds on 
the elliptieity of the equator as sensed by Syncom II in a seven month drift period over Brazil to 
be 


J 22 = - (1.9 1.2) x 10-® 


(17a) 


or 


a 0 - b Q = 73 + 8 meters , 

which is the difference in major and minor equatorial radii of the earth (Reference 1), and 


\ 22 = - (20 ± 7 ) ° with respect to Greenwich , (17b) 

which locates the major axis of the elliptical equator of the earth. 
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CONCLUSIONS 


From the evidence of many recent geoids based on different kinds of gravity data the following 
conclusions are drawn: 

1. Earth gravity reductions of long term drift accelerations of 24-hour satellites can, at 
favorable longitudes, ignore all but the second-order effect without incurring large errors. 

2. Estimates of earth equatorial ellipticity based on Syncom n drift over Brazil sensitive to 
only the second order effect were: 

J 22 = - (1.70 ±0.05) x io' 6 : 

\ 22 = - (19+6)° • 

When one accounts for the probable range of all higher order gravity effects on the 24-hour 
satellite, on the basis of lower altitude data only, these figures should be 

J 22 = - (1.9 ±0.2) x 10’ 6 

and 

X 22 = - ( 20 ± 7)° . 

3. Further data from Syncom n should clarify the relatively small uncertainty left in exact 
knowledge of the earth's elliptical equator. 

4. A very conservative estimate of the absolute maximum tangential velocity requirements 
for station keeping of 24-hour equatorial satellites is (see Appendix C) 

AV t = 7.5 ft./sec./yr. 


(Manuscript received June 3, 1965) 
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Appendix A 


The Longitude Gravity Field of the Earth 
According to Investigations Since 1915 


Table A-l lists 30 sets of longitude coefficients for the earth’s gravity potential as reported by 
geodesists since 1915. The significance of the gravity constants with respect to the shape of the 
sea-level geoid is discussed in Reference 1 and illustrated in Reference 8. 
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Longitude Coefficients in the Earth’s Gravity Potential Jv E - £ £ £ [l - (R 0 /r )" P n m ( s i n J nm cos m (\ - \ nm )] l , as Reported 1915-19642. 



16 



References 


"Determination of the Ellipticity of the Earth’s Equator from 
Observations on the Drift of the Satellite" NASA Technical Note 
D-2759, May 1965. 

Private communication from W. M. Kaula (July 1964). 

Private communication from W. M. Kaula (July 1964). 

Private communication from W. M. Kaula (July 1964). 

Guier, W. H., and Newton, R. R., "Non-Zonal Harmonic Coefficients 
of the Geopotential from Satellite Doppler Data," Silver Spring, Mary- 
land: Johns Hopkins University Applied Physics Laboratory, TG-520, 
November 1963. 

Kaula, W. M., "Improved Geodetic Results from Camera Observations 
of Satellites," J. Geophys. Res. 68(18):5183-5190, September 15, 1963. 

Izsak, I. G., "Tesseral Harmonics in the Geopotential," Nature, 
199(4889):137-139, July 13, 1963. 

Private communication from W. M. Kaula: (May 1963). 

Kaula, W. M., "Tesseral Harmonics of the Gravitational Field and 
Geodetic Datum Shifts Derived from Camera Observations of Satel- 
lites," J. Geophys. Res. 68(2):473-484, January 15, 1963. 

Private communication from W. M. Kaula (July 1964). 

Private communication from Y. Kozai: (Oct. 1962). 

Private communication from Newton: (April 1962). 

Newton, R. R., "Ellipticity of the Equator Deduced from the Motion of 
Transit 4A," J. Geophys. Res. 67(1):415-416, January 1962. 

Kozai, Y., "Tesseral Harmonics of the Gravitational Potential of the 
Earth as Derived from Satellite Motions,” Astron. J. 66(7):355-358, 
September 1961. 

Kaula, W. M., "A Geoid and World Geodetic System Based on a Com- 
bination of Gravimetric, Astrogeodetic and Satellite Data," J. Geophys. 
Res. 66(6):1799-1811, June 1961. 

Izsak, I. G., "A Determination of the EUipticity of the Earth's Equator 
from the Motion of Two Satellites,” Astron. J. 66(5):226-229, June 1961. 


(19) Kaula, W. M., "Analysis of Satellite Observations for Longitudinal 

Variations of the Gravitational Field," In: Space Research II, ed. by H. C. 
van de Hulst, C. de Jager, and A. F. Moore: 360-372, Amsterdam: 
North-Holland Publishing Co., 1961. 

(20) Nicolaides, J. D. and Macomber, M. M., "Passive Dynamics in Space 

Flight," Bureau of Naval Weapons Paper, 1962. 

(21) Kaula, W. M., "Statistical and Harmonic Analysis of Gravity," J. 

Geophys. Res. 64(12):2401-2421, December 1959. 

(22) Jeffreys, H., The Earth: its Origin, History, and Physical Constitution, 

4th ed., London: Cambridge University Press, 1959:127-194. 

(23) , (26), (28), (29), (30) Heiskanen, W. A., and Vening-Meinesz, F. A., The Earth and its 

Gravity Field, New York: McGraw-Hill Book Co., 1958:79. 

(24) Zhongolovitch, I. D., "Potential of the Terrestrial Attraction," Akad. 

Nauk SSSR Inst. Theoret. Astron. 6(8):505-523, 1957 (In Russian). 

(25) Subbotin, M. F., Kurs Nebesnoi Mekhaniki, Vol. 3, 2nd ed., Moscow: 

Gostekhizdat, 1949. 

(27) In: Monthly Notices Royal Astronom. Soc. Geophys. Suppl., 5(55):1942. 


18 



Appendix B 


Calculation of Apparent Equatorial Eliipticity From Syncom II - Orbit 
Drift Over Brazil, in the Gravity Field of Kaula - Combined (1964) 

From Appendix A the gravity constants of the geoid of Kaula - Combined (1964) are 


J 22 

= “ 1.51 x 

VO 

1 

0 

r*H 

^33 

- 22.8° 

^22 

= - 15.5° 


J 42 

= -0.163 x 10" 6 

J31 

~ “l.Slx 

10'® 

^42 

= 37.0° 

^31 

= 0.0° 


J 44 

= -0. 0053 x 10' 

J 33 

= -0. 149 X 

10'® 

^44 

= 35.8° . 


Thus, from these values in Equation 9a, if Syncom II had been drifting in the gravity field of this 
geoid, at k = - 54.8° its acceleration would have been 

1(- 54.8°) Kaula . Coiib 19M - 5. 06 (- 1. 51 x 10'®) sin 2( - 54. 8° + IS. 5°) 

-0.0605 (~ 1. 51 x 10'®) sin (- 54.8° -0.0°) + 5. 26 (-0.149x 10'®) sin 3(- 54.8°- 22.8°) 

-0.0158 (-0.163x 10'®) sin 2(- 54- 8° -37.0°) +6.83 (-0.0053x 10'®) sin4(- 54.8° -35.8°) 

= 6.793 x 10'* , (B-l) 

in units of - 12 n 2 (R 0 /a s ) 2 rad/day 2 . Similarly, at -59.2°, the acceleration would have been 

M-59.2°) Kaula . Coaib W64 = 5.06 (- 1. 51 x 10'® ) sin 2( _ 59. 2° + 15. 5°) 

-0.0605 (- l.Slx 10'®) sin (- 59.2° -0.0°) + 5.26 (-0. 149 x 10"®) sin 3(- 59.2°- 22.8°) 

-0.0158 (-0. 163 x 10~« ) sin 2(- 59.2° -37.0°) +6.83 (-0.0053x 10'®) sin 4(- 59.2°- 35.8°) 

= 6.815x10 ® , (B-2) 

in units of - 12 n 2 (R 0 /a s ) 2 rad/day 2 . 
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Under the assumption that these accelerations arise from the second order field only, (the 
assumption made in Reference 4) from Equations 13a, B-l, and B-2, the apparent equatorial major 
axis location is calculated at 


f - (-0.87965 X 6. 793 + 0.94206 X 6.815) 1 
*22 2 tan ( (-0.33545x 6.815+0.47562x 6.793) J 


r -12.6° . (B-3) 

It is seen that the error in \ 22 made under the simplifying assumption for this geoid is 

A ^ 22 = \ J2 (actual) - A. ' (apparent) 


= - 15.5 -(- 12.6) = -2.9° (B-4) 

(see Table 1). Under the same simplifying assumption, from Equations 13b, B-l, and B-2, an ap- 
parent magnitude of equatorial ellipticity is calculated from Syncom n drift above this geoid, as 


- 1 . 292 x 10“ 6 


(-0. 87965 x 6. 793+0. 94206 x 6.815) 2 
+ (-0. 33545 x 6. 815+0. 47562 x 6.793) 2 


1/2 


= - 1.35 x 10“ 6 . 


(B-5) 


Expressed as a ratio, the error in j 22 made under the simplifying assumption for Syncom n drift 
above this geoid, is 


lu _ - 1.51x 10~ 6 

j 22 - 1 . 35 X 10" 6 


1.12 . 


(B-6) 


See Table 1 for the results of this simulated drift calculation for seven other recent geoids based 
on different kinds of gravity data (as listed in Table A-l). 
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Appendix C 


Probable Maximum East-West Station Keeping Velocity 
for 24-Hour Equatorial Satellites 


On the basis of an equatorial ellipticity with a gravity constant given by 

J„ = “1.7 x 10“ 6 , (C-l) 

it was calculated in Reference 4 that the maximum tangential station keeping requirement for 24- 
hour equatorial satellites for positions midway between the major and minor equatorial axes (i.e., 
at about 26°, 116°, -154° and -64°) is 


AV t - 5.36 ft/sec-yr , (C-2) 

which is proportional to J 22 . This result was predicted on the basis of an entirely negligible in- 
fluence of higher than second-order earth longitude gravity on 24-hour satellites. From the section 
of this report covering the probable range of higher order distortions to equatorial ellipticity the 
most conservative estimate of the probable distortion in the above Brazil - Syncom II j 22 determin- 
ation due to higher order effects is seen to be 

J 22 = - (1. 9+0.2) X 1(T 6 = -2.1 xlO' 6 . (C-3) 

R is pointed out in Reference 7 that absolute maximum accelerations on 24-hour satellites in a 
fourth order field are likely to be as much as 12% higher at certain longitudes than the maximum 
accelerations in a simple second order field. 

A conservative 15% margin on the velocity from the above maximum value of j 22 gives a 
very conservative value of 


AV t (max. ) = * 15 * * 5- _ 7 5 ft/sec-yr (C-4) 

that may be necessary for station keeping of 24-hour satellites. From Figure 1 it is seen that the 
longitude location for this absolute maximum requirement is most likely to be over Indonesia in 
the vicinity at 116° East. 
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Appendix D 

List of Symbols 


a 0 ,b 0 Major and minor radii of the elliptical equator of the earth (according to the convention 
in Reference 1, R 0 =(a 0 + b 0 )/2) 

a s Semi major axis of the synchronous, or 24-hour satellite 

F( i) nm Inclination factor for the nm gravity component expressing the reduction in mean long 
term longitude acceleration of an inclined 24-hour satellite over that for an equatorial 
24-hour satellite 

Jnm’ Gravity constants of the earth’s field; magnitude and phase angle of the rm spherical 
gravity harmonic (Reference 1) 

J 22 * ^22 Apparent values of the second-order earth longitude gravity constants computed in the 
Syncom n - Brazil experiment under the simplifying assumption that the higher order 
longitude gravity field is negligible (See text) 

r 0 Mean equatorial radius of the earth (- 6378.2 km) 

r 22 , Z ^ 22 Bias factors of j 22 and k 22 resulting from the simplifying assumption that the higher 
order longitude gravity field is negligible (See text) 

AV X Tangential velocity requirement for perfect station keeping of an equatorial "synchronous" 
satellite 

K Geographic longitude east of Greenwich 

Earth’s Gaussian gravity constant (~ 3.986 X 10 5 km 3 /sec 2 ) 


NASA- Langley, 1966 
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